Ground-state phase diagram of the two-dimensional t-J model 
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The ground-state phase diagram of the two-dimensional t-J model is investigated in the context of the tensor 
network algorithm in terms of the graded Projected Entangled-Pair State representation of the ground-state wave 
functions. There is a line of phase separation between the Heisenberg anti-ferromagnetic state without hole and 
a hole-rich state. For both J = OAt and J = 0.8r, a systematic computation is performed to identify all the 
competing ground states for various dopings. It is found that, besides a possible Nagaoka's ferromagnetic state, 
the homogeneous regime consists of four different phases: one phase with charge and spin density wave order 
coexisting with a p,.(/? y )-wave superconducting state, one phase with the symmetry mixing of d + s-wave 
superconductivity in the spin-singlet channel and p. v (p v )— wave superconductivity in the spin-triplet channel in 
the presence of an anti-ferromagnetic background, one superconducting phase with extended .s-wave symmetry, 
and one superconducting phase with p. v (p v )-wave symmetry in a ferromagnetic background. 
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Since the discovery of high temperature superconductiv- 
ity id, significant progress has been made in our understand- 
ing of strong correlation physics. It was Anderson [2] who re- 
alized the importance of Mott-Hubbard insulators and put for- 
ward the resonating valence bond picture as a promising route 
towards the understanding of an electron pairing mechanism 
responsible for an unprecedented high transition temperature, 
which are observed for copper oxides (cuprates). Actually, 
the detailed analysis of electronic states deduced from exper- 
iments shows that the undoped parent compound is a Mott- 
Hubbard insulator and the hole doping is mainly on oxygen 
sites, with its effective low energy physics described by the 
two-dimensional t-J model ^ . 

Although a lot of efforts have been made to gain a full pic- 
ture of the underlyi ng p hysics of the two-dimensional t-J 
model (see, e.g., B4I 417IO . no consensus has been achieved as 
to the question whether or not the two-dimensional t-J model 
superconducts. On the one hand, the variational Monte Carlo 
(VMC) method clearly indicates that the d^t_y2~wave super- 
conductivity is stable at absolute zero temperature for a phys- 
ically realistic coupling parameter regime 11811 . and field the- 
oretical slave-boson approximation yields qualitatively many 
peculiar phenomenological features of cuprate superconduc- 
tors fl9tl . On the other hand, exact diagonalization (ED) of 
the t-J model on a small cluster and quantum Monte Carlo 
(QMC) simulation of the two-dimensional Hubbard model do 
not produce convincing evidence supporting the existence of 
superconductivity in the t-J model within a physically rele- 
vant regime |8l l20fj . 

In this paper, we attempt to investigate the ground-state 
phase diagram of the two-dimensional t-J model by exploit- 
ing a newly-developed tensor network algorithm 12111 in terms 
of the graded Projected Entangled-Pair State (gPEPS) repre- 
sentation of the ground-state wave functions (for an ungraded 
version, see 12211 ). It is found that, the algorithm yields reli- 
able results for the two-dimensional t-J model at and away 
from half filling, with the truncation dimension up to 6. We 
are able to locate a line of phase separation (PS) between the 



Heisenberg anti-ferromagnetic state without hole and a hole- 
rich state, which qualitatively agrees with the results based on 
the high-temperature expansions (HTE) [11], the VMC lfl3[l. 
and the density-matrix renormalization method (DMRG) 11411 . 
In the homogeneous regime, the two-dimensional t-J model 
exhibits very rich physics. Away from half filling, the regime 
is divided into four different phases: one phase with charge 
and spin density wave order coexisting with a p. v (p v )— wave 
superconducting state (CDW+SDW+PW), one phase with 
the symmetry mixing of d + s-wave superconductivity in 
the spin-singlet channel and p x (p y )—v?aye superconductiv- 
ity in the spin-triplet channel in the presence of an anti- 
ferromagnetic background (DSW+PW+AF), one supercon- 
ducting phase with extended s-wave symmetry (SW), and 
one superconducting phase with p x (Py)— wave symmetry in a 
ferromagnetic background (PW+FM), besides a possible Na- 
gaoka's ferromagnetic state (FM) 12311 . 

Ground-state phase diagram of the two-dimensional t-J 
model. The two-dimensional t-J model is described by the 
Hamiltonian 18[] consisting of a hopping term and a super- 
exchange interaction: 

H = -t J] VKcicjo- + H.c.)P] + J £(S; • Sj - \tn nj), (1) 



(0> 



<«7> 



where S; are spin 1/2 operators at site i on a square lattice, P is 
the projection operator excluding double occupancy, and t and 
J are, respectively, the hoping constant and anti-ferromagnetic 
coupling between the nearest neighbor sites (ij). 

The model is simulated by exploiting the gPEPS tensor net- 
work algorithm [21], with the truncation dimension up to 6. 
Here, we stress that, although the algorithm is variational in 
nature, an assumption has to be made on the choice of the 
unit cell of the gPEPS representation of the ground-state wave 
functions. For the t-J model, we have chosen the plaquette 
as the unit cell (see Fig.Q] left panel). 

In Fig. |2] we plot the ground-state phase diagram for the 
two-dimensional t-J model, with the vertical and horizon- 
tal axes, n and J It, denoting the number of electrons per site 
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FIG. 1: (color online) Left panel: The unit cell of the graded Pro- 
jected Entangled-Pair State representation of the ground-state wave 
functions. Right panel: The pattern of the coexisting charge and spin 
density wave order. Here, the radius of the circles are proportional 
to the fillings at sites, whereas the arrows inside the circles represent 
the directions and magnitudes of the spin density wave order param- 
eters. This represents a vertical commensurate stripe state, which 
breaks the four-fold rotation symmetry and the translation symmetry 
in the horizontal direction for charge density wave order and in both 
directions for spin density wave order. 



and the ratio between the anti-ferromagnetic coupling and the 
hoping constant, respectively. At half filling, n = 1, the t - J 
model reduces to the two-dimensional Heisenberg spin 1 /2 
model. Therefore, a long range anti-ferromagnetic order ex- 
ists iU In this case, our algorithm yields the ground-state 
energy per site, e = -1.16757, for the truncation dimension 
D = 4, and e — -1.16837, for the truncation dimension 
D = 6. This is comparable to the best QMC simulation re- 
sult: e — -1.16947 [24, 25], although the anti-ferromagnetic 
Neel order moment is 0.36 versus 0.31. Away from half fill- 
ing, the model exhibits quite different behaviors for small 
and large values of the anti-ferromagnetic coupling 7. For 
7/f > 0.95, there is a line of PS between the Heisenberg 
anti-ferromagnetic state without hole and a hole-rich state, 
whereas for 7/f < 0.95, no PS occurs. This agrees qual- 
itatively with the results based on the HTE ifllll . the VMC 
method [13], and the DMRG [14]. Note that our result for 
the transition point J c = 3.45? at low electron density is quite 
close to the exact value 7 C = 3. 4367 f lfl5ll . Here, the simula- 
tion has been performed for the truncation dimension B = 6. 

We point out that a discrepancy exists concerning the PS 
boundary of the t - 7 model. In Ref. 0], a combination of 
analytic and numerical calculations is used to establish the ex- 
istence of PS at all super-exchange interaction strengths. Al- 
though the simulation based on the Green's function Monte 
Carlo method supports this scenario fioll . many others JTII - 
[3 found that the model phase separates only for 7/f larger 
than some finite critical value around 7/f ~ 1. That is, PS oc- 
curs only outside the physically realistic parameter region of 
the f - 7 model. With the observation that no significant shift 
with the truncation dimension B increasing from B = 4 lElll 
to B = 6 is found, we conclude that PS does not occur for 
J It < 0.95. 

In the homogeneous regime, the two-dimensional f - 7 
model exhibits very rich physics. Away from half filling, both 
the d + s-wave superconductivity in the spin-singlet channel 




FIG. 2: (color online) The proposed ground-state phase diagram of 
the two-dimensional t — J model. First, for J/t > 0.95, there is a 
line of phase separation (PS), whereas for J/t > 0.95, no PS oc- 
curs. Second, the homogeneous regime is divided into four different 
phases: one phase with charge and spin density wave order coexist- 
ing with a p x (p v )-wave superconducting state (CDW+SDW+PW), 
one phase with the symmetry mixing of d + s-wave superconduc- 
tivity in spin-singlet channel and p_ v (p v )-wave superconductivity in 
spin-triplet channel in the presence of an anti-ferromagnetic back- 
ground (DSW+PW+AF), one superconducting phase with an ex- 
tended s-wave symmetry (SW), and one superconducting phase with 
a p-wave symmetry (PW+F), besides a possible Nagaoka's ferro- 
magnetic state (FM). Here, a systematic computation has been per- 
formed for both J/t = 0.4 and J/t = 0.8 with the truncation dimen- 
sion D = 6, whereas the dash lines separating different phases are a 
guide for the eyes. For J/t = 0.4, the DSW+PW+AF phase occurs 
for fillings from n = 0.818 (denoted as point a) to n = 0.765 (de- 
noted as point b), the SW phase occurs for fillings from n = 0.765 
to n = 0.087 (denoted as point c), and the CDW+SDW+PW and 
PW+FM phases occur for fillings from n = 1 to n = 0.818 and 
from n = 0.087 to n = 0, respectively. For J/t = 0.8, the 
DSW+PW+AF phase occurs for fillings from n = 0.877 (denoted 
as point d) to n = 0.665 (denoted as point e), the SW phase occurs 
for fillings from n = 0.665 to n = 0.051 (denoted as point f), and the 
CDW+SDW+PW and PW+FM phases occur for fillings from n = 1 
to n = 0.877 and from n = 0.051 to n = 0, respectively. 



and the p x (py)— wave superconductivity in the spin-triplet 
channel, and/or charge and spin density wave order, occur 
in different doping regimes for a fixed 7/f. Here, a super- 
conducting state is characterized by a superconducting order 
parameter A = (A), with A defined as follows: For s-wave, 
A s = l/(4V2) ( -; ,.[<•:.„ -<■: ,,..+<•;.,.,.-<•:., 

]; for ./-wave, A d = 1/(4 VI) ^[c^+c+^-c^- 

ct . _ u ] - [f<->4]; for p x -wave, A Pl = A Pjt+ - A Pl -, with A Px± = 

1/2 (ct . T c+ , . T ,[c + . T c+ . . +c + . ,c+ .J/ V2,c + . ,ct . . X 

and a similar definition of A p for p y -wa\e. 

Charge and spin density wave order coexisting with a 
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Px(Py)— wave superconducting state. For a physically realistic 
super-exchange coupling J and a hopping constant t (such as 
J It — 0.4), a spin-triplet p v (p v )-wave superconducting state 
coexists with charge and spin density wave order for dopings 
up to around 6 ~ 0.18, with 5 = 1 - n. The occurrence of 
such a spin-triplet p v (p v )-wave superconducting component 
is unexpected, although this would, in our opinion, have been 
anticipated from the presence of the Nagaoka's ferromagnetic 
state 12311 . In addition, Dagotto and Riera 18J observed im- 
portant ferromagnetic and anti-ferromagnetic correlations in 
this regime from the ED of the t - J model on a small cluster, 
which may be properly interpreted as the precursor of a spin- 
triplet p x (p v )-wave pairing state, coexisting with charge and 
spin density wave order. Note that the coexisting charge and 
spin density wave order forms a pattern, as displayed in Fig.Q] 
(right panel). In this phase, the translational invariance under 
one-site shifts, and four-fold rotation symmetry, as well as the 
t/(l) symmetry in the charge sector and the S U(2) symmetry 
in the spin sector, are spontaneously broken. 

In Fig. [3] the magnitude of the spin-triplet p x (p y )-wave su- 
perconducting order parameter A p , along with those of the 
charge and spin density wave order parameters, are plotted 
for both J/t — 0.4 and Jjt — 0.8, which have been evaluated 
with the truncation dimension D = 6. In this phase, all the 
order parameters are evaluated for two sites A and B in the 
unit cell, with the subscripts A and B as their labels, except 
for the charge density wave order parameter that is defined as 
half the difference between the densities at two sites A and B 
(see Fig. [H right panel). 

The mixing of the spin-singlet d + s—wave and spin-triplet 
Px(Pv)—wave superconductivity in the presence of an anti- 
ferromagnetic background. As shown in Fig. 12 there is a 
superconducting phase with mixed spin-singlet d + s-wave 
and spin-triplet p x (p y )-wa\e symmetries in the presence of 
an anti-ferromagnetic background. For / = 0.4, it occurs for 
fillings from n = 0.818 to n — 0.765. For J = 0.8, it occurs for 
fillings from n = 0.877 to n — 0.665. In this phase, the trans- 
lational invariance under one-site shifts, and four-fold rotation 
symmetry, as well as the t/(l) symmetry in the charge sector 
and the S U (2) symmetry in the spin sector, are spontaneously 
broken. Note that the of-wave, s-wave and p x (p y )-v?axe com- 
ponents are homogeneous. 

In Fig. [3] the magnitudes of the cZ-wave, s-wave and 
Px(Py)— wave order parameters, Aj, A.„ and A p (upper panel), 
along with that of the anti-ferromagnetic order parameter 
(lower panel), are plotted for the t—J model with J/t = 0.4 and 
J/t = 0.8 which have been evaluated with the truncation di- 
mension D = 6. One observes that the of-wave component Ad 
vanishes at electron fillings n c \\ n c \ = 0.765 and n c i = 0.818 
for J/t = 0.4 and n cl = 0.665 and n c2 = 0.877 for J/t = 0.8. 
We also list the numerical values of the d-wave and .s-wave 
components, Ad and A s , for both J/t = 0.4 and J/t = 0.8 at 
different fillings in Table U together with their ratio A s /A,/. 
The fact that the pairing symmetry (in the spin-singlet chan- 
nel) is of the d + s-wave nature manifests itself in that the ratio 
A s /Ad is always real. On the other hand, the symmetry mix- 



J/t 


n 


A., 


A,, 


AJA d 


0.4 


0.768 


0.0107 - 0.0042/ 


-0.0256 + 0.0100/ 


-0.416 


0.776 


0.0101 - 0.0040; 


-0.0261 + 0.0103/ 


-0.386 


0.785 


0.0096 - 0.0037/ 


-0.0266 + 0.0103/ 


-0.360 


0.794 


0.0091 -0.0032/ 


-0.0270 + 0.0100/ 


-0.336 


0.804 


0.0085 - 0.0033/ 


-0.0270 + 0.0105/ 


-0.314 


0.814 


0.0080-0.0031/ 


-0.0272 + 0.0106/ 


-0.295 


0.8 


0.720 


0.0147 - 0.0052/ 


-0.0375 + 0.0134/ 


-0.391 


0.722 


0.0145 - 0.0053/ 


-0.0378 + 0.0133/ 


-0.383 


0.725 


0.0143 - 0.0050/ 


-0.0382 + 0.0133/ 


-0.375 


0.741 


0.0129 - 0.0045/ 


-0.0398 + 0.0139/ 


-0.325 


0.757 


0.0116-0.0046/ 


-0.0401 + 0.0160/ 


-0.289 


0.800 


0.0091 -0.0037/ 


-0.0419 + 0.0168/ 


-0.217 


0.811 


0.0085 - 0.0034/ 


-0.0421 + 0.0169/ 


-0.203 


0.833 


0.0076 - 0.0030/ 


-0.0422 + 0.0169/ 


-0.179 


0.842 


0.0072 - 0.0029/ 


-0.0419 + 0.0168/ 


-0.171 



TABLE I: The numerical values of the s-wave and rf-wave supercon- 
ducting order parameters, A s and Ad, and their ratio A s /Aj for both 
J/t = 0.4 and J/t = 0.8 at different fillings in the DSW+PW+AF 
phase. 



ing of the spin-singlet and spin-triplet channels arises from the 
spin-rotation symmetry breaking: spin is not a good quantum 
number. 



The spin-singlet superconducting phase with extended 
s—wave symmetry. For a fixed 0.1 < J/t < 0.95, if one keeps 
on increasing doping until the cZ-wave component Aj van- 
ishes, then only an s-wave order parameter A s survives, re- 
sulting in a spin-singlet superconducting phase with extended 
s-wave symmetry. For J/t = 0.4 and J/t — 0.8, the s-wave 
pairing order parameter A s is plotted in Fig. [3] A peculiar fea- 
ture of the s-wave order parameter is that it is almost linearly 
increasing with increasing doping, until it reaches its maxi- 
mum, and then it monotonically decreases. 

Although this phase is perhaps only of academic interest 
due to its unphysical nature of large dopings for cuprate su- 
perconductors, we believe it is important to clarify its pairing 
symmetry. Note that a small portion located between 2.0 < 
J/t < 3.4367 has been identified as a spin-singlet supercon- 



ducting phase with s-wave symmetry in Refs. (91 15, 26ll . 

77ie spin-triplet superconducting phase with p x (p y )—wave 
symmetry in a ferromagnetic background. For large dopings, 
a spin-triplet superconducting phase with p x (p y )-wave sym- 
metry in a ferromagnetic background occurs. The existence of 
such a spin-triplet superconducting phase with p x (p y )-wa\e 
symmetry has been pointed out by Kagan and Rice l26ll . 
Therefore, our simulation results are consistent with their an- 
alytical analysis. 

In Fig. [3] the magnitude of the spin-triplet p-wave order 
parameter, together with that of the ferromagnetic order pa- 
rameter, are plotted for both J/t = 0.4 and J/t = 0.8 with the 
truncation dimension D = 6. Here, the p x (p y )-wa\e order 
parameter is identical at all sites, thus it is homogeneous. 

Doping-induced quantum phase transitions. Now we turn 
to doping-induced quantum phase transitions (QPTs) for the 
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FIG. 3: (color online) Upper panel: The dependence of the rf-wave, 
.s-wave, and p v (p v )-wave superconducting order parameters, A d , A s 
and A,,, on electron filling n. Here, we evaluate A s , Aj and A p for 
J/t = 0.4 and J/t = 0.8. Lower panel: The dependence of the spin 
density wave (SDW) and the charge density wave (CDW) order pa- 
rameters in the CDW+SDW+PW phase, the anti-ferromagnetic Neel 
order parameter in the DSW+PW+AF phase, and the ferromagnetic 
order parameter in the PW+F phase on electron filling n. Note that, 
in the CDW+SDW+PW phase, all the order parameters are evaluated 
for two sites A and B in the unit cell, with the subscripts A and B as 
their labels, except for the charge density wave order parameter that 
is defined as half the difference between the densities at two sites A 
and B (see Fig.[T] right panel). 



two-dimensional t — J model. For J/t = 0.4, the model 
undergoes a QPT from the CDW+SDW+PW phase to the 
DSW+PW+AF phase at n = 0.818, and a QPT from the 
DSW+PW+AF phase to the SW phase at n = 0.765, and 
again a QPT from the SW phase to the PW+F phase at 
n = 0.087. For J/t = 0.8, the model undergoes a QPT from 
the CDW+SDW+PW phase to the DSW+PW+AF phase at 
n = 0.877, and a QPT from the DSW+PW+AF phase to the 
SW phase at n — 0.665, and again a QPT from the SW phase 
to the PW+F phase at n = 0.051. In addition, a QPT oc- 
curs between the DSW+PW+AF phase and the PW+F phase. 
Here, all QPTs are first-order, as one may judge from the be- 
haviors of the order parameters shown in Fig. [3] In passing, 
we point out that the above results have been deduced from 
the ground-state fidelity approach in the context of tensor net- 
work representations 1 27], with superconducting order param- 



eters being read off from both one-site and two-site reduced 
density matrices according to a general scheme advocated in 
Ref. Hi . 

In summary, we have investigated the ground-state phase 
diagram of the two-dimensional t - J model in the context of 
the tensor network algorithm. The relevance of our simulation 
results to the high temperature cuprate superconductors will 
be discussed in a forthcoming publication M29I1 . 

This work is supported in part by the National Natural 
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